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QCD resummation for jet substructures 
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We provide a novel development in jet physics by predicting the energy profiles of light-quark 
and gluon jets in the framework of perturbative QCD. Resumming large logarithmic contributions 
to all orders in the coupling constant, our predictions are shown to agree well with Tevatron CDF 
and Large-Hadron-Collider CMS data. We also extend our resummation formalism to the invariant 
mass distributions of light-quark and gluon jets produced in hadron collisions. The predicted peak 
positions and heights in jet mass distributions are consistent with CDF data within uncertainties 
induced by parton distribution functions. 
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PACS numbers: 12.38.Cy,12.38.Qk,13.87.Ce 

It has been a long-standing challenge to predict sub- 
structures (including energy profiles and masses) of light- 
quark and gluon jets in the perturbative QCD (pQCD) 
theory. During the Tevatron run 1 era in the early 1990s, 
it was found that next-to-leading-order (NLO) QCD cal- 
culations cannot describe experimental data on jet sub- 
structures. Hence, it has been a custom for experimental- 
ists to compare jet substructures measured at the Teva- 
tron, either at run 1 or run 2, with predictions from full 
event generators such as PYTHIA or HERWIG. While 
the full event generators (usually with specific tuning) 
could describe data, it remains desirable to develop a 
theoretical framework for the study of jet substructures. 
For that, the soft-collinear effective theory was adopted 
in the literature, such as [l|,|2j. In this Letter, we propose 
a novel approach to predicting jet substructures based on 
the pQCD resummation formalism [3j. We show that re- 
sults of the resummation formalism for light-quark and 
gluon jets are well consistent with the energy profiles 
measured by CDF at Tevatron |4j and CMS at Large 
Hadron Collider (LHC) |5[, and with the mass distribu- 
tions measured by CDF [6(. 

It is known that the top quark is predominantly pro- 
duced at rest at the Tevatron and can be clearly iden- 
tified by detecting three (or more) isolated jets from its 
hadronic decays. However, this strategy will not work for 
identifying a highly boosted top quark [7H10l| at the LHC, 
which results in a single jet. Furthermore, it has been 
pointed out llj, |l2| that an energetic QCD (light-quark 
or gluon) jet can have an invariant mass around the top- 
quark mass to fake a top-quark jet. That is, a boosted 
top-quark jet is difficult to be discriminated against an 
ordinary QCD jet. This difficulty also appears in the 
identification of a highly boosted Higgs boson decaying 
into a bottom-quark pair [Lj, [lj] , or a highly boosted 
W or Z boson decaying into hadronic final states, for 
they can all produce a single-jet experimental signature. 
In order to improve jet identification at the LHC, ad- 
ditional information from jet energy profiles is needed, 
because jets initiated by different parent particles usu- 



ally produce different energy profiles. 



We denote 
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the jet energy function as 
Jf(Mj, Pt,v 2 , R, r) for defining a light-quark (/ = q) or 
gluon (/ = g) jet with mass Mj, transverse momentum 
Pt, and cone size R, which describes the all-order energy 
distribution within a smaller cone of size r < R. J? 
is constructed by inserting a sum of the step functions 
Si kiT<d(r — 9i) into the usual jet definition [15j, where 
kiT and Oi are the transverse momentum and the angle 
of the final-state particle i with respect to the jet 
axis. At leading-order (LO), it is a S- function, i.e., 
Jt — PtS(Mj), which is independent of r, because 
6 = at this order. The jet definition contains a 
Wilson line along the light cone, which collects gluons 
emitted from other parts of the collision process and 
collimatcd to the parent particle of the jet. To employ 
the resummation technique, wc vary the Wilson line 
into an arbitrary direction n^ with n 2 ^ [16j |. The 
dependence of Jf on n^ and the jet momentum Pj 



Mj, and Pj ■ n 



appears through the invariants n 2 , P 2 
(which is related to Pt). When r approaches zero, the 
phase space of real radiation is strongly constrained, 
so the infrared enhancement in real radiation does not 
cancel completely with that in virtual correction. The 
resultant large logarithms of the ratio (Pj ■ n) 2 /(n 2 r 2 ), 
which is conveniently defined as \P? P T j (Ar 2 )\v 2 , should 
be resummed to all orders in the coupling constant a B . 
It is easy to see from the above ratio that the variation in 
n, i.e., v 2 , can turn into the variation in r. To compare 
with present experimental data on jet energy profile , we 
consider the jet energy function with the jet invariant 
mass being integrated out, which corresponds to taking 
the N = 1 moment in the Mellin space and is denoted 
as J?(1,Pt,v 2 ,R,r). By definition, different choices of 
n^ yield the same collinear divergences associated with 
the jet. Hence, the effect of varying n, i.e., v 2 , does not 
involve the collinear divergences and can be factorized 
out of the jet energy function, leading to an evolution 



equation 
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The one-loop kernel G^ (K^) absorbs the hard (soft) 
dynamics of the variational effect, whose expressions are 
similar to those derived in [17J , but with the step function 
in angle being inserted into the real-gluon piece of K^ 1 ' . 
The solution describing the evolution of the jet energy 
function from the initial value uf n = C\r 2 /(C 2 R 2 ) to the 
final value v\ = 1 is written as 

jf(l,Pr,4,R,r) = Jf(l,P T ,vLR,r) 
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The color factor Cf is equal to Cf(= 4/3) and Ca(= 3) 
for the light-quark and gluon jet, respectively, /3q is the 



QCD Beta function [18(, and n/ is the number of active 
light-quark flavors. The value of i/f n diminishes the large 
logarithms in the initial condition J e (1,Pt, tf n , R,r), 
which is then evaluated up to NLO including non- 



logarithmic-r terms. The value 



1 implies the pres- 



ence of the large logarithms in Jf(l, Pt, v$, R, r), which 
have been summed into the Sudakov integral in Eq. ©. 
We set the 0(1) constants C\ = Ci = 1 and C = 
cxp(5/2) (C = exp(17/6)) for the quark (gluon) jet in or- 
der to reproduce the large logarithms a s In r and a s In r 
in the NLO calculations. The variation of these 0(1) 
constants reflects theoretical uncertainty in our formal- 
ism. The value of r in the lower bound is taken to be 
larger than 0.1, so that it is safe to evaluate the Sudakov 
integral perturbatively. We then derive the energy pro- 
file ^l(f) |4{ as the energy fraction accumulated within 
the cone of size r < R in terms of the solution in Eq. © , 
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which respects the normalization ^(r = R) = 1. 
Equation ((4]) contains the convolution of the LO differ- 
ential cross section daf /dPr and the quark and gluon jet 
energy functions. Using the CTEQ6L parton distribu- 
tion functions (PDFs) [19j , we compare the resummation 
and NLO predictions in Fig. Q] with the Tevatron CDF 
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FIG. 1: Resummation (solid) and NLO (dashed) predictions 
for jet energy profiles compared with CDF data [J]. 



data [J]. The agreement between the resummation pre- 
dictions and the CDF data is obvious for all Pt values. 
As Pt increases, the accumulation of energy inside the 
jets becomes faster. The NLO predictions derived from 
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FIG. 2: Resummation (solid) and NLO (dashed) predictions 
for jet energy profiles compared with CMS data Q]. 



Jf (1, Pt, v\, R, r) are also displayed, which overshoot 
the data. Figure [5] shows the agreement of the resum- 
mation predictions for the jet energy profiles with the 
LHC CMS data at 7 TeV @| and the overshooting of the 
NLO predictions. The above consistency indicates that 
our resummation formalism has captured the dominant 
dynamics in a jet energy profile and can give a direct and 
reliable prediction for this observable. 

As stated before, the invariant mass distribution of an 
energetic jet can be utilized as an experimental signa- 
ture of new physics. We shall demonstrate below that 
our formalism is also applicable to this jet substruc- 



ture by predicting the mass distributions of the light- 
quark and gluon jets. Following the previous analysis 
on jet energy profiles, we vary the Wilson line into an 
arbitrary direction ?i M with n 2 ^ 0, when implement- 
ing the resummation technique [16| to derive the evo- 
lution equation for the jet function Jf(Mj,PT,v 2 ,R). 
In this case, the overlap of the collinear and soft en- 
hancements generate the double logarithms of the ratio 
(Pj ■ n) 2 /(M 2 jn 2 ) = {R 2 P T /(AMj))u 2 , so that the varia- 
tion of n, i.e., v 2 , can turn into the variation of Mj. The 
dependence on the jet cone size R is introduced through 
the Mellin transformation with respect to the variable 
x = Mj/(R 2 P^). The solution describing the evolution 
of the jet function from the initial value v 2 n = Ci/^C^TV), 
where N = N exp(7^) with 7^ being the Euler's con- 



stant, to the final value 



1 is given by 



term (as suggested by S(N) in Eq. ((6])), as well as a linear 
term in N (for describing a final-state jet [20[). The pa- 
rameters ao,i.2 can be determined by a fit to experimen- 
tal data for certain jet momentum and jet cone. In this 
work, we perform fits to full event generators PYTHIA 
211 ] and SpartyJet [22j for the quark (gluon) jet produced 
at the Tevatron run 2 energy 1.96 TeV with P T = 600 
GeV and R = 0.7, which yields a = -0.35, a x = 0.50 
(—4.59), and ai = —1.66. With this nonperturbative 
contribution, we are ready to predict jet mass distribu- 
tion for arbitrary values of Pr and R using the improved 
resummation solution: 



Jf*>(N,Pr,i4,R) 



J f (N,P T ,v 2 n ,R) 
xexp[S(N) + S NP {N)] (10) 



J f (N,P T ,vlR) = J f (N,P T ,v? n ,R)exv[S(N)], (5) 



The jet function in Mj space is derived via the inverse 
Mellin transformation 



in which the Sudakov exponent is written as 

S(N) = - f 2 ^L\ A { as (y 2 R 2 P 2 ))\J^ 
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The value v\ = 1 implies the presence of the large log- 
arithms (in powers of In TV) in J/(iV, Pr,i/|,i?), which 
have been summed into the Sudakov exponent. The 
initial condition Jf(N, Pr, v 2 ^ R) can then be evaluated 
perturbatively. At NLO, they have the expressions 
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for the light-quark and gluon jets, respectively. The re- 
maining TV-dependent terms are suppressed by 1/N and 
their effect is small. 

For a fixed Pt, the scale RPt/N involved in Eq. © 
becomes so low at extremely large N that the pertur- 
bative analysis fails and nonperturbative contributions, 
arising from hadronization and underlying events, need 
to be included in order to predict jet mass distribution 
in the small Mj region. For convenience, we introduce 



S» r (N) 



N 2 



R 2 P 2 



j{C f a \nN + a 1 ) + C f a 2 



NQg 
RP T ' 
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into the Sudakov exponent with Qo being set to 1 GeV. 
It consists of a logarithmic term and a Gaussian smearing 



1 



-N 7RES/ 



J^(Mj) = — / dN(l - x)-»Jf**(N), (11) 
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where x = M 2 /(R 2 P T ) and the contour C runs along 
the negative real axis in the complex- N plane and circles 
the origin counterclockwise with a finite radius. To avoid 
the Landau pole, we flatten the running behavior of a s 
18|,l23| at certain scale /i c : 



a s O) 



a s (^ c exp[iArg(/i)]), |/x| < /x c 
Oi t (n), \fi\ > He 



(12) 



In the small Mj region, all moments in N are equally im- 
portant, and those containing powers of In TV dominate. 
Therefore, our resummation formalism, in which these 
large logarithms are summed up, can give reliable predic- 
tions. In the large Mj region, the nonlogarithmic terms 
are not negligible, so we have to improve the resumma- 
tion formula by including the nonlogarithmic terms (say, 
up to NLO), namely, by matching fixed-order results to 
the resummation formula as done for the Drell-Yan pro- 
cesses [24j. 

The jet mass distribution is calculated by convoluting 
the LO differential cross section of inclusive (quark or 
gluon) jet production with the corresponding quark or 
gluon jet function, 



da 
dMj 



f 



dP T ^fj f (M 2 : P T ,4,R). (13) 



We are now ready to compare our predictions for the jet 
mass distributions with the Tevatron and LHC data, by 
choosing C\ = exp(7s), C2 = exp(— je) and \x c = 0.3 
GeV in the numerical analysis with the CTEQ6L PDFs 
[19( . Results from different choices of C± , C'2 and \x c re- 
veal theoretical uncertainty, which will be investigated 
in a forthcoming paper. The comparison to the Teva- 
tron data [6j, with the kinematic cuts Pt > 400 GeV 



and 0.1 < \Y\ < 0.7, is presented in Fig. [3j where the 
label NLL/NLO denotes the prediction of NLL resum- 
mation with the NLO initial conditions given in Eqs. (O 
and IJHJ). It shows that our predictions agree well with 
the CDF data for R = 0.4 and 0.7 in the intermediate 
jet mass region. In the small invariant mass region, the 
predictions can describe the overall shape of the jet dis- 
tributions within the uncertainty induced by PDFs |6j, 
though the peak positions are slightly lower than the 
CDF data. However, for large M.j, e.g. Mj > 100(200) 
GeV for R = 0.4(0.7), the resummation prediction drops 
off quickly and deviates from data. This can be further 
improved by matching to exact calculations of NLO and 
beyond. The resummation predictions for the jet mass 



CDF 6 fb' 1 , R=0.4 
NLL/NLO, R=0.4 
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FIG. 3: Comparison between resummation predictions and 
Tevatron data. 



and jets from boosted Higgs boson and weak gauge bo- 
son (W, Z or Z', etc.) hadronic decays will be studied 
in a forthcoming paper. We shall also calculate other 
jet substructures, such as the angularity [1|. These stud- 
ies are crucial for the LHC physics program in terms of 
testing the QCD theory and identifying new physics sig- 
nals. We have also applied the resummation formalism to 
predict the light-quark and gluon jet mass distributions. 
To describe jet mass distributions in the low mass region, 
we need to introduce some nonperturbative contributions 
in the resummation formalism. The relevant parameters 
can be fixed at one jet energy scale, and then employed 
to make predictions for other energy scales, which are 
found to agree with the Tevatron data. It is expected 
that our formalism will successfully apply to upcoming 
LHC data of jet mass distributions. 
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distributions at Tevatron with R = 0.3 and at LHC with 
R = 0.7 are shown in Fig. @J which can be tested by 
future measurements. 




FIG. 4: Resummation predictions for jet mass distributions 
at Tevatron and LHC. 

In conclusion, we have developed a theoretical frame- 
work based on the pQCD theory for analyzing the sub- 
structures of the light-quark and gluon jets. This is the 
first time in the literature that pQCD is shown to de- 
scribe well the jet energy profiles and mass distributions, 
which are the most commonly discussed physical observ- 
ables to describe the substructure of a jet signal at hadron 
colliders. We have demonstrated that the resummation 
predictions for the jet energy profiles, with the jet invari- 
ant mass being integrated out, are insensitive to nonper- 
turbative inputs, and in excellent agreement with both 
the Tevatron CDF and LHC CMS data, for arbitrary 
values of jet momentum. In view of the fact that the jet 
energy profile is a useful feature for jet identification at 
the LHC, the energy profiles of boosted top-quark jets, 
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